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Abstract. 

A detailed analysis of the radiative and Auger de-excitation channels of K-shell vacancy states in Fe n-Fe ix has been carried out. 
Level energies, wavelengths, A-values, Auger rates and fluorescence yields have been calculated for the lowest fine-structure 
levels populated by photoionization of the ground state of the parent ion. Different branching ratios, namely Ka 2 P^a t , Kfi/Ka, 
KLM/KLL, KMM/KLL, and the total K-shell fluorescence yields, a> K , obtained in the present work have been compared with 
other theoretical data and solid-state measurements, finding good general agreement with the latter. The Ko^/Ko-i ratio is found 
to be sensitive to the excitation mechanism. From these comparisons it has been possible to estimate an accuracy of ~10% for 
the present transition probabilities. 
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1. Introduction 

The iron K lines appear in a relatively unconfused spectral re- 
gion and have a high diagnostic potential. The study of these 
lines has been encouraged by the quality spectra emerging from 
Chandra and by the higher resolution expected from Astro-E 
and Constellation-X. In addition there is a shortage of accurate 
and complete level-to-level atomic data sets for the K-vacancy 
states of the Fe isonuclear sequence, in particular for the lowly 
ionized species. This undermines line identification and real- 
istic spectral modeling. We are currently remedying this situa- 
tion by systematic calculations using suites of codes developed 
in the field of computational atomic physics. Publicly avail- 
able packages have been chosen rather than in-house develop- 
ments. In this context, complete data sets for the n — 2 K- 
vacancy states of the fi rst row, namely Fe x vrn -Fe xxv, have 
been reported earlier by Bau tista et alJ ll20C)3h and lPalmeri et alJ 
to be referred to hereafter as Paper I and Pa per II, and 
for the second row (Fe x-Fe xvn) bv lMendoza et alJ J2003I) . to 
be referred to as Paper III. 

The K lines from Fe species with el ectron occupancies 

N > 17 have been studied very little. Ijacobs & Rozsnval 

(1986) have computed fluorescence probabilities in a frozen- 
core approximation for vacancies among the n£ subshells of 
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the Fe isonuclear sequence. Kaastr a"& Mewd i 19931) have cal- 
culated the inner-shell decay processes for all the ions from Be 
to Zn by scaling published Auger and radiative rates in neu- 
trals. Both of these studies ignore multiplets and fine-structure. 
Otherwise, the bulk of the data in the literature is devoted to 
solid-state iron. Regarding wavelengt hs and line intensit y ra- 
tios, numerous references are listed in Hol zer et al J (1 19971) who 
measure the Kai 2 and emission lines of the 3d transi- 

tion metals using a high-precision single-crystal spectrometer. 
Fewer p ublications are availa ble on the experimental K Auger 
spectra: iKovalik et alJ dl987b have derived the KLM/KLL and 
KMM/KLL ratios from the Auger electron spectr um of iron 
produced by 57 Co decay, and lNemethv et all i 19961) have mea- 
sured the KLL and KLM spectra of the 3d transition metals but 
have not determined the KLM/KLL ratio. Concerning K-shell 
fluorescence yields, measuremen ts covering th e period 1978— 
93 have been reviewed by Hubbell et al. ( 1994) following ma- 
jor compilations by Bambvn ek et alJ(ll972l) and lKrausd (11979). 

The present work is a detailed analysis of the radiative and 
Auger de-excitation channels of the K-shell vacancy states in 
the third-row species Fe 11-Fe ix. Energy levels, wavelengths, 
A-values, Auger rates and fluorescence yields have been com- 
puted for the lowest fine-structure levels in configurations ob- 
tained by removing a Is electron from the ground configura- 
tion of the parent ion. In Section 2 the numerical method is 
briefly described. Section 3 outlines the decay trees and se- 
lection rules. Results and discussions are given in Section 4, 
followed by a summary and conclusions (Section 5). All the 
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atomic data calculated in this work are available in the elec- 
tronic Tables 3—4. 

2. Numerical method 

In Paper III we conclude that the importance of core-relaxation 
effects increases with electron occupancy, so calculations of 
third-row iron ions require a computational platform which is 
well suited to treating these effects. For this reason the cal- 
culatio ns reported he re are carried out using the hfr pack- 
age b v lCowanl dl98ll) although auto structure jBadnelll ll986, 
1997) is also heavily used for comparison purposes. In hfr 
an orbital basis is obtained for each electronic configuration 
by solving the Hartree-Fock equations for the spherically av- 
eraged atom. The equations are derived from the application 
of the variational principle to the configuration average en- 
ergy and include relativistic corrections, namely the Blume- 
Watson spin-orbit, mass-velocity and the one-body Darwin 
terms. The eigenvalues and eigenstates thus obtained are used 
to compute the wavelength and A-value for each possible tran- 
sition. Autoionization rates are calculated in a perturbation the- 
ory scheme where the radial functions of the initial and final 
states are optimized separately, and configuration interaction 
is accounted for only in the autoionizing state. Configuration 
interaction is taken into account among the following config- 
urations: (3d + 4s) M , [3p](3d + 4s) M+1 , [2p](3d + 4s) M+1 and 
[Is] (3d + 4s) M+1 where [nl] stands for a hole in the nl sub- 
shell and (nl + n'l') M represents all possible distributions of M 
electrons among the nl and n'l' shells, M ranging from M = 
in Fe ix to M — 7 in Fe n. 

Given the complexity of [2p] [3p] double-vacancy channels 
in ions with an open 3d shell, the level-to-level computation of 
the Auger rates with hfr and autostructure proved to be in- 
tractable. However, average values can be used for all the lev- 
els taking advantage of the near constancy of the total Auger 
widths in ions for which the KLL channels reduce the outer 
configuration to a spectat or (see Paper III). The refore, we em- 
ploy the formula given in Palm eri etall feoOl) for the single- 
configuration average (SCA) Auger decay rate 



,SCA 



n"l"sl"') 



(1) 



Zi(2Jj + 1) 

where the sum runs over all the levels of the autoionizing con- 
figur ation, s/s is a statistical factor given in Eqs. (15-16) of 
Pal meri et alJ J200 H) that contains the dependence on the active 
shell (nl, n'l', n"l") occupancy, wdA a (nln'l' -> n"l"sl"') is 
the two-electron autoionization rate which is a function of the 
radial integral s and for which the c omplete expression is given 
in Eq. (11) o flPalmeri et al.l |2001 ). The calculated Auger rates 
are expected be as accurate as those obtained in a level-to-level 
single-configuration approach. 

3. Photoionization selection rules and decay trees 

We have focused our calculations on the Fe K-vacancy states 
populated by photoionization of the ground state of the parent 
ion 



where // is the outer configuration of the ground state, i.e. 
3d 6 4s 2 in Fe i, 3d 6 4s in Fe n and 3d N (N - 1 - 6) in Fe m- 
Fe vm, [Is] denoting a K hole. In Eq. Q, K-shell phot oion- 
ization leads to a p wave with selection rules specified bv lRavJ 
Jl976h : U = L, S' - S = ±1/2 and J' - J = +1/2, +3/2. 
Therefore, only few states are expected to be populated. 

The radiative and Auger decay manifolds of a K-vacancy 
configuration [Is]// can be outlined as follows: 



- Radiative channels 
[Is]// — > [3p]yu + y/j 

— > [2p]/z+7a 

- Auger channels 
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(2S + 1) 



Lj + y — -> [is]// 



(2S' + 1) 



L' r +, 
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where the negative exponent of // stands for the number of elec- 
trons that have been extracted from the outer subshells. In the 
radiative channels, forbidden and two-electron transitions have 
been excluded as it has been confirmed by calculation that 
they display very small transition probabilities (logA r < 11). 
Therefore, the two main photo-decay pathways are the result 
of the 2p — > Is and 3p — » Is single electron jumps that give 
rise respectively to the Ka (~ Al. 94) and K/3 (~ Al. 75) arrays. 

It has also been numerically verified that the participator 
Auger channels (KMM(p) and KLM(p)) contribute less than 
one percent to the total Auger widths, and hence, they have not 
been taken into account. Of primary interest in the present work 
is the branching ratios of the Ka, K/3, KMM, KLM and KLL 
channels and their variations with electrons occupancy N. 

4. Results and discussion 

Centroid wavelengths for the Ka and K/3 unresolved transition 
arrays (UTAs) computed with hfr are presented in Tabled in- 
cluding also a comparison with experiment. For the Kai and 
Ka2 lines, it can be seen that the agre ement with the solid-state 
measurements bv lHolzer et alJ ( fl9 97) (~ 0.5 mA ) is somewhat 
better than that with the EBIT results for Fe x JPecaux et alJ 
1995) of within 2 mA, and the slight blueshift with increas- 
ing predicted by hfr is consistent with experiment. A small 
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Fig. 1. Comparison of centroid wavelengths for (a) Ka and (b) 
KjS UTAs in Fe ions with 17 < N < 25. Fil led circles: this 
work. Filled triangles: Kaastra & Mewel (ll993l) . 



redshift with N is also found for the KJ3 array. On the other 
hand, as shown in Figure ^ the present findings contrast with 
the steeper redsh i ft for both UTAs obtained from the data by 
lKaastra&Mewe1 ( ll993l) . 

In Table|2] Ka^/Ka^ intensity ratios are tabulated for differ- 
ent ionization stages. Two cases have been considered: HFR1, 
only the decay lines from the K-vacancy levels populated by 
the photoionization of the respective ground state are included; 
and HFR2, all the transitions from the levels belonging to the 
[ls](3d + 4s) M+1 complex are taken into account. Noticeable 
differences in the ratios from these two cases may be ap- 
preciated indicating sensitivity to the excitation mechanism. 
Th e HFR2 r atio fo r Fe n is closer to the Dirac-Fock value 
of IScofieldl (1974), who also considered all the K-vacancy 
states, and to the solid-state experiments (fk olzer e{ alj|l997 



iMcGrarv et alJl97lUSalem & Wimm eJ 1 97(AIWilliam sil 5 
In Figure |2] theoretical and experimental Kfi/Ka intensity 
ratios are plott ed as a function of electron number. Most experi- 
mental ratios (Rao et al. 1986; Berenyi et al. 1978; Salem et al. 
Il972l ISlivinskv & Ebertlll972l lHansen et alJll97(j) have been 
scaled down by 8.8% in order to extract the radiative-Auger 
(5%) and K/? 5 (3.8%) satellite co ntributions fro m the KyS 
UTA which theory does not include ( Verma 2000). The value 



Fig. 2. Comparison of Kfi/Ka intensity ratios for Fe ions 
with 18 < N < 25. Filled circles: hfr, this work. 
Dotted upright triangles: auto structure, this work. Dotted 
diamonds : mcdf- s al, th is work. Filled upright triangles: 
lKaastra & Mewel dl993h . Filled square s: Lfacobs & Rozsnvail 
dl986l). Filled diamonds: IScofieldl i 1974 . Filled inverted trian - 
gles: IJankowski & Polasikl(ll989|) . Circle s: IPeruio et al.Nl987l) 
Uprigh t triangles: iHolzer et alJ Jl997h. Squares: iRao et all 
dl986l). Diamonds : iBer^nvTet alJ dl978l). Inverted tr i angles : 




Fig. 3. Comparison of Auger rates for Fe ions with 18 < 
N < 25: (a) KLM/KLL and (b) KMM/KL L ratios. Filled cir - 
cles: hfr , this work. Upright f il led tri angles: IChen e t alJ dl97^>. 
Circle s: iBhalla &~R amsdale lll97oj). Squares: [ Kovalik et alJ 
Jl987h . Upright triangles: iMehlhorm & Albridgellll963l) 
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Table 1. Comparison of centroid wavelengths (A) for the Ka 
and K/5 unresolved transition arrays in Fe ions with 17 < N < 
25. The hfr wavelengths have been weighted with the fluores- 
cence yields. Experime ntal w avelengths for Fe x (N = 17) 
are from iDecaux et alJ d 19951) and for solid-state iron from 
lHolzeretal.ldl997l) . 



HFR 



Expt 



N 


Ka 


KB 


17 


1.9369 


1.7488 




1.9407 


1.7493 


18 


1.9366 


1.7511 




1.9402 




19 


1.9364 


1.7527 




1.9403 




20 


1.9362 


1.7540 




1.9405 




21 


1.9360 


1.7557 




1.9404 




22 


1.9357 


1.7566 




1.9400 




23 


1.9355 


1.7565 




1.9397 




24 


1.9356 


1.7565 




1.9398 




25 


1.9356 


1.7566 




1.9398 





Ka 



KB 



1.9388(5) 
1.9413(5) 



1.936041(3) 
1.939973(3) 



1.756604(4) 



Table 2. Comparison of hfr Ka^/Kai ratios for Fe ions (18 < 
N < 25) with experiment and previous theoretical estimates. 



N 


Ion 


HFR1" 


HFR2 


18 


Fe ix 


0.7000 


0.5912 


19 


Fe vm 


0.5780 


0.6015 


20 


Fe vii 


0.4541 


0.5821 


21 


Fe vi 


0.5545 


0.5538 


22 


Fe v 


0.5852 


0.5346 


23 


Fe iv 


0.4426 


0.5143 


24 


Fe m 


0.4438 


0.5032 


25 


Fe ii 


0.4476 


0.4969 



Expt. 



Theory 



0.51(2) c 
0.507(10) d 
0.506 6 
0.4998 f 



0.5107s 



a Includes only the decay lines from the K-vacancy levels listed in 
Table 3 

b Incl udes all the decay lines from the [ls](3d + 4s) M+1 complex 
c iHolzer et alJll997h 



McGrarv et all i 197 II) 
Salem & WimmeiHl970h 
WmamsNl933l) 



Scofiel 



nsR1933D 
d 119741) 



quoted in IPeruio et all i ll 987b does not take into account the 
blend with the Kfis satellite and has been corrected accord- 
ingly. With the exception of the values by iKaastra & Mewd 
(1993), theory predicts a decrease of the ratio with N, and 
the theoretical scatter is comparable with that among the solid- 
state experiments of just under 20%. The computed results by 
hfr, autostructure and Lfacobs & Rozsnvaj {l986) for N = 




Fig. 4. Comparison of K-shell fluorescence yields, u>k, for Fe 
ions 10 < N < 25. Filled circles: hfr value obtained for the 
first K-vacancy level, 10 < N < 17 (Paper III), 18 < N < 25 
(this work). Filled squares: autostructure value obtained for 
the first K-vacancy level, Paper III. Filled upright triangle: 
hfr value obtained by averaging on all the K-vacancy levels, 
10 < N < 17 (Paper III), 18 < N < 25 (this work). Filled in- 
verted triangle: autostructure value obtained by averagi ng on 
all the K-vacancy levels, Paper III. Circles: Oacobs & Rozsnvail 
dl986l) . Squares: IKaastra & Mewd dl993l). Filled diamonds: 
value recommended by Hubbell et al. ( 1994). Upright triangle: 
I Sold d 19921) Inve rted triangle: IPious et al.1 dl992l) . Diamond: 
Bhan et all dl98 lb . Dotted circle: Sing h et alJ 



25 are in g ood agreement w i th the bulk of the experimen- 
7 ■ II — I II 1 I 

tal values (Pernio et al. 1987; Rao et al. 1986; Berenvi et al 



1978; Salem et al. 1972; Sli vinskv & Eberil972tlHansen et al 
1970). Holzer et al. ( 1997) mention a possible systematic de- 
viation in one of their corrections to explain the discrepancy 
of their data with ot her measuremen ts. On the o ther han d, the 
Dirac- Fock ratios by Scofieldl d 19741) and l.Tankowski & Polasik 
(1989) at N = 25 appear significantly higher. We have there- 
fore proceeded to veri fy these values by us ing the same code 
(mcdf-sal) as in IJankowski & Polasikl |l989), and as shown in 
Figure 13 they are accurately reproduced; moreover, a similar 
decrease in the ratio with N is also obtained with this method. 
The spread of the different data sets in this comparison indi- 
cates a probable accuracy of our hfr transition probabilities of 
~ 10%. 

A decrease with N is also predicted by hfr for the Auger 
KLM/KLL and KMM/KLL ratios (see Figure For Fe n, 
pr esent results are in ex cellent accord with the measurements 
of iKovalik et alJ d 19871) and other theory dChen et all fr 9791: 
Bhalla & Ramsdal dl 19701) but somewhat higher than the older 
experimental estimate of Mehlhorm & Albridge ( 1963). The 
total K-shell fluorescence yields, u>k, for Fe ions with 10 < 
N < 25 are presented in Figure|4] In both hfr and autostruc- 
ture data sets, the fluorescence yields have been computed for 
fine-structure K-vacancy levels: 



u> K (i) = 



A r {i) 



A r (i) + A a (i) 



(10) 



In order to display the effect of the level population, we report 
in Figure[3]the value for the first K-vacancy level and the yield 
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Fig. 5. Simulations carried out near 6.4 keV with the xstar modeling code for several values of the ionization parameter 



averaged over all the K-vacancy levels. As expected, the flu- 
orescence yield is independent of the population mechanism 
in the third row ions where the Kfi, KLM and KMM decay 
channels become spectators. Given the complexity of level-to- 
level calculations of Auger rates, it has not been possible to 
compute fluorescence yields with autostructure for N > 17. 
The most recent ex perimental measurem ents in the s olid ( Sole 
ll992HPious et alll992tlsingh et alll990tlBhan et alll98lh are 
plotted along with the value recommended by iHubbell et all 
( 1994). They show a scatter within 10% and are slightly lower 
than our hfr value f or N - 5, i.e. 5 to 1 0% with respect to 
IHubbell et all (1 1994 and lBhan et ail dl98ll) respectively. This 
is not a significant discord considering the expected accuracy 
of our fluorescence yields (~ 20%) as discussed in Paper III. 
As mentioned in the previous section, the missing participator 
Auger channels will not contribute significantly to the Auger 
widths (less than 1%). Concerning the shake-up channels, they 
will not affect the Auger widths because they are due to mix- 
ing between valence shell configurations which have the same 
Auger widths. We have verified that this is the case in a sim- 
ple calculation ([ls1(3d + 4s) 8 -> [2p] 2 (3d + 4s) 8 +/fcs). Data by 
iKaastra & Mewd i ll 9931) have been corrected so as to reproduce 



the best solid-state yields compiled in Bambvn ek et alJ i 19721) 
and are therefore somewh at lower than our values. Results by 
Ijacobs & RozsnvaH i 19861) for N > 12 predict a substantially 
lower yield. 

In order to evaluate the impact of the new atomic data on 
Fe K-line modeling, the data sets generated in the present work 
and those in Paper s I— III have been inc luded in the xstar mod- 
eling code jKallman & Bautistall200ll) . Runs for different ion- 
ization parameters, are shown in Figure|5] The other param- 
eters have been assigned the following values: cosmic abun- 
dances; a gas column density of 10 16 crrT 2 ; a gas density of 
10 12 cirT 3 ; an X-ray source luminosity of 10 38 erg/s; a power- 
law index of 1 ; and energy bins (or channel widths) of 1 eV. In 
Figure [5] one clearly sees that the UTA centroid is redshifted 
when £ goes from 0.001 to 1 and then blueshifted for £ > 1 . The 
shape of the UTA changes also considerably for ^ > 1 where 
the ionization balance favors the first and second row ions. 

Electronic Tables 3—4 list radiative and Auger widths for 
295 energy levels; and wavelengths, A-values and fluorescence 
yields for 396 transitions. 
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5. Summary and conclusions 

Following the findings of our previous study o n the se cond-row 
iron ions (Paper III), the hfr package (Cowan 1981) has been 
selected to compute level energies, wavelengths, decay rates 
and fluorescence yields for the K-vacancy states in Fe n-Fe ix. 
The calculations have focused states populated by ground- 
state photoionization. Due to the complexity of the level-to- 
level Auger cal culations, we have employed a compact formula 
( Pal meri et all200ll) to compute Auger widths from the hfr ra- 
dial integrals. 

The hfr centroid wavelengths for the Ka a nd KB UTAs 
in Fe n reproduce the solid-state measurements (IHolzer et alJ 
1 19971) to better than 1 mA. Moreover, the red shift pre- 
dicted by hfr for the Ka lines in species with higher ioniza- 
tion stage is in acc ord with the EBIT wavelengths in Fe x 
JPecaux e t al.|l99 5|) thus contra dicting the previous trend spec- 
ified bv lKaastra & Mewd dl993l) . We have also carried out ex- 
tensive comparisons of different hfr branching ratios, namely 
Ka 2 /Ka u K/3/Ka, KLM/KLL, KMM/KLL and cj k , with other 
theoretical and experimental data. The present ratios for Fe n 
are in good agreement with the solid-state measurements, and 
the Kq'2/Kq'i ratio has been found to be sensitive to the excita- 
tion mechanism. It has been possible from these comparisons 
to estimate an accuracy of ~10% for the hfr transition proba- 
bilities. 

The new atomic data sets for the whole Fe isonuclear se- 
quence that have emerged from the present project are provid- 
ing a more reliable platform for the modeling of Fe K lines. 
Preliminary si mulations of the emissivit y of a photoionized 
gas with xstar ( Kallman & Bautista 2001) have shown Ka line 
profiles and wavelength shifts that are sensitive to the ioniza- 
tion level of the gas. Further work will therefore be concerned 
with improving diagnostic capabilities by means of accurate K- 
shell photoionization and electron impact excitation cross sec- 
tions. 
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